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R. N.  A.  Plimmer 


SUMMARY 

A  theoretical  analysis  is  presented  of  the  number  of  launches  that  are 
required  to  establish  and  maintain  a  satellite  communication  system.  A  simple 
model  is  used  to  represent  the  process  in  that  launches  are  assumed  to  take 
place  with  a  certain  a  priori  probability  whilst  the  satellites  have  an 
exponential  probability  distribution  of  lifetimes. 

Graphs  and  equations  are  given  for  the  various  probability  distributions 
of  the  system  during  the  establishment  and  maintenance  phases  of  the  system  and 
the  effect  of  the  time  interval  between  launches  is  also  examined  by  using  a 
Markov  chain  representation  of  the  system.  It  is  noted  that  many  of  the 
probability  distributions  are  well  represented  by  a  Gaussian  distribution. 

The  results  quite  clearly  indicate  the  manner  in  which  the  number  of 
launches,  and  hence  the  cost,  depend  critically  upon  the  launch  success 
probability  and  the  ratio  of  the  launch  interval  to  the  satellite  mean 
lifetime. 
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1  INTRODUCTION 

Several  civil  and  military  satellite  communication  systems  have  been 
proposed  in  recent  years  whereby  radio  signals  for  telecommunication  purposes 
are  transmitted  between  ground  stations  via  orbiting  satellites.  In  order  to 
provide  a  communication  channel  which  is  available  at  any  time  of  "016  day  these 
systems  require  several  satellites  to  be  continually  maintained  in  orbit ; 
satellites  which  have  failed,  as  regards  communication  capabilities,  must  be 
replaced  by  further  satellite  launchings  so  as  to  maintain  sufficient  satellites 
in  orbit.  The  number  of  satellites  required  in  orbit  may  vaiy  from  a  minimum 
of  three  for  a  geostationary  system  to  between  fifty  and  a  hundred  for  a  random 
system.  Since  the  running  costs  for  a  random  system  are  high,  principally  on 
account  of  the  large  number  of  satellites  required,  such  a  system  seems  rather 
unattractive  at  the  moment.  Consequently,  as  far  as  this  Report  is  concerned, 
investigation  has  been  restricted  to  medium  and  high  altitude  station  keeping 
systems  where  not  more  than  about  20  satellites  are  required  in  orbit. 

Of  the  satellites  launched,  only  a  certain  proportion  will  be 
operationally  satisfactory,  in  that  they  will  be  injected  successfully  into  the 
correct  orbital  position  with  all  their  mechanical  and  electrical  systems 
working.  Present  day  estimates  suggest  that  the  probability  of  such  a 
successful  launch  may  be  as  low  as  3Q/o  although  no  doubt  this  will  be  improved 
with  future  experience  and  research  on  the  reliability  of  components.  Even  when 
a  fully  operating  satellite  has  been  launched  successfully  into  orbit,  it  will 
eventually  fail  to  function  due  to  deterioration  or  breakdown  of  its  numerous 
working  components  and  will  have  to  be  replaced  by  another  satellite.  Present 
design  of  satellite  components  aims  at  achieving  a  mean  lifetime  of  about  five 
years  and  eventually  ten  years. 

For  costing  purposes,  preplanning  of  launch  vehicle  pads,  components,  and 
facilities,  it  is  essential  to  know  the  likely  number  of  satellite  launchings 
and  the  possible  launching  rates  which  will  be  required  in  order  to  establish 
and  maintain  the  system  in  operation.  In  this  Report  an  attempt  has  been  made 
to  assess  these  quantities  on  the  basis  of  a  relatively  simple  model  of  the 
system.  It  is  also  shown  that  many  of  the  probability  distributions  may  be 
well  approximated,  for  practical  purposes,  by  Gaussian  distributions,  thus 
enabling  rapid  estimates  of  system  performance  to  be  made  without  detailed 
calculations. 
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2  MATHEMATICAL  MODEL  OF  THE  SYSTEM 

Having  decided  upon  a  particular  satellite  communication  system,  let  it 
be  assumed  that  a  total  of  at  least  N  orbiting  satellites  are  required  in  order 
to  provide  continuous  communications.  Thus  N  satellites  must  first  be  established 
in  their  respective  orbits  and  then,  when  any  satellite  fails,  it  must  be  replaced 
as  quickly  as  possible  by  another.  Since  a  satellite  cannot  be  replaced  immediately, 
it  is  evident  that  more  than  N  satellites  must  be  maintained  in  orbit  to  obtain  a 
high  probability  of  continuous  communication  or  else,  perhaps  for  economic  reasons, 
an  appreciable  break  in  the  communication  link  must  be  tolerated.  One  possible 
method  of  attempting  to  achieve  a  continuous  communication  system  is  simply  to 
duplicate  the  satellites  in  pairs  and  switch  to  the  remaining  satellite  when  its 
neighbour  fails.  Although  a  slight  break  in  transmission  is  still  inevitable  at 
switching  it  will  not  mean  a  complete  breakdown  unless  the  neighbour  has  also 
failed.  Even  so,  such  a  method  would  only  seem  feasible  for  a  stationary  system, 
for  example,  attempting  to  maintain  six  satellites  (three  pairs)  in  orbit.  For 
medium  altitude  systems  demanding  a  greater  number  of  satellites,  a  more  sensible 
compromise  would  be  to  maintain  one  or  two  spare  satellites  in  the  system  and 
manoeuvre  them  into  position  following  the  failure  of  a  satellite  and  accepting 
the  inevitable  break  in  communication.  The  communication  time  lost  by  this 
method  would  most  likely  be  less  than  in  attempting  to  replace  the  lost  satellite 
by  a  new  launching.  Rerouting  of  the  communication  link  following  a  satellite 
failure  may  also  be  possible  if  the  increased  transmission  delays  could  be 
tolerated. 

In  order  to  formulate  a  mathematical  model  for  representing  the  mechanism 
of  establishing  and  maintaining  a  communication  system  certain  simplifying 
assumptions  have  been  made  here.  In  the  first  place,  no  account  has  been  taken 
of  the  particular  configuration  of  the  communication  system  considered,  i. e. 
whether  polar,  equatorial,  stationary,  constant  ground- track  etc,  so  that  a 
communication  system,  for  the  purposes  of  this  Report,  simply  consists  of  a 
number  of  satellites  N  which  must  be  established  and  maintained  in  orbit. 

Secondly,  as  regards  the  launching  of  satellites,  the  location  of  the  launching 
pads  and  time  taken  to  place  a  satellite  into  its  correct  station  do  not  enter 
into  the  analysis.  In  fact  the  launching  system  is  considered  to  be  capable 
simply  of  launching  satellites  into  orbit  at  a  certain  rate  and  thus  the  number 
of  available  launching  pads  only  determines  the  frequency  with  which  satellites 
may  be  launched.  Furthermore,  no  account  is  taken  of  the  variation  of  firing 
rates  between  different  launching  pads  and  the  probability  of  successfully 
launching  an  operationally  satisfactory  satellite  into  its  correct  position  in 
orbit  is  assumed  to  have  a  constant  value  denoted  by  p. 
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Finally  since  the  operational  lifetime  of  a  satellite  is  determined  by 
the  correct  functioning  of  numerous  pieces  of  equipment  it  seems  very  reasonable 
to  assume  that  the  lifetime  of  a  satellite  has  a  Poisson  probability  distribution* 
Thus  we  will  assume  that  the  probability  that  a  satellite  fails  in  the  small 
increment  of  time  dt  between  times  t  and  t  +  dt  is  dt/£,  independent  of  the  time 
t  and  where  £  is  the  mean  lifetime  of  the  satellite*  Consequently,  the 
probability  that  the  satellite  operates  for  a  time  t  and  then  fails  in  the 
following  interval  dt  is 


p(t)  dt  s  e  (l ) 

and  the  probability  that  it  functions  for  at  least  time  t  is 

oc 

P(t)  =  f  &  =  e-V*  .  (2) 

t 

It  may  be  noted  in  passing  that  this  type  of  probability  distribution  is 
particularly  convenient  mathematically  since  the  initial  operational  time  of  a 
satellite  does  not  enter  into  the  analysis  simply  because  the  probability  of 
any  functioning  satellite  failing  in  the  following  instant  is  independent  of 
time.  Any  other  probability  distribution  would  greatly  complicate  the  analysis 
since  the  distribution  would  depend  upon  the  launch  time  of  the  satellite. 

3  ESTABLISHMENT  OF  A  SATELLITE  SYSTEM  USING-  HI&H  LAUNCHIN&  RATES 

3.1  Theoretical  probability  distribution 

Let  us  assume  that  the  probability  of  failure  of  ary  satellite  in  the 
system  is  negligible  over  the  established  phase,  or,  in  other  words ,  that  the 
launching  rate  is  extremely  high.  In  this  limiting  case,  the  number  of 
satellites  which  have  to  be  launched  depends  solely  upon  the  probability  of 
making  a  successful  launch,  p.  Now,  the  probability  p®(n,p)  that  n  firings 
will  have  to  be  made  before  there  are  N  successful  launches  is  equal  to  the 
probability  that  the  last  firing  is  successful  times  the  probability  of  (N-1 ) 
successes  in  the  previous  (n-1 )  firings.  Hence 

Pfl(njp)  =  P  (n-1,  p)  (3) 

where  (n,p)  is  the  binomial  distribution  function 

Bn  (n,p)  =  f  „)  p  (1-p) 
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and  thus 

P^(n;p)  =  (-)pVp)”-N  (5) 

is  the  required  probability  density.  Consequently,  the  probability  of  obtaining 
N  satellites  in  orbit  using  no  more  than  n  launchings  is  simply 


p>;p> 


pVp) 


r-N 


r=N 


(N  *  1)'  (6) 


since  each  term  of  this  series  represents  the  probability  of  the  event  that  the 
Nth  success  occurs  on  the  rth  firing  and  all  these  events  (r  =  N,  N  +  1,  ...  n) 
are  mutually  exclusive.  Equation  (6)  thus  gives  the  probability  distribution 
function  of  the  number  of  satellites  successfully  launched  into  orbit  as  a 
function  of  the  number  of  firings. 

An  alternative  expression  for  P.jfT(n;p)  can  be  found  from  the  following 

JM 

reasoning:  in  n  firings  the  probability  of  N  successful  launchings  is  B^(n,p) 
and,  consequently,  in  n  firings  the  probability  of  at  least  N  successes  is 


n 


n 


P^(n;p) 


Br(n,p) 


Y  (r) 


r=N 


(7) 


The  equivalence  of  equations  (6)  and  (7),  straightforward  to  prove  algebraically, 
is  evident  from  the  reasoning  leading  to  these  expressions  since  either  contains 
all  the  possible  ways  of  obtaining  N  satellites  in  orbit  without  making  more 
than  n  attempted  launches. 


3.2  Evaluation  of  the  probability  distribution 

As  they  stand,  the  expressions  (6)  and  (7)  are  in  an  extremely  inconvenient 
form  for  calculation  purposes  but  they  may  be  transformed  into  alternative 
tabulated  forms.  The  Taylor  expansion  of  a  function  f(x)  may  be  written  with 
remainder  in  the  form 


f(q+p) 


(q+tp)  dt 


(8) 


0 
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and  thus,  taking  f  r=  (q+p)*1  with  q  =  1-p,  this  becomes 


N-1 


1  = 


1 

n\  r, .  \n-r  f  n!  N/ ,  ,  \N-1  r* 

r)  p  +  J  (k-1  ) !  (n-N)!  p  ^  [l  ‘ 


p(l-t)]n"K  dt 


r=0 


Changing  the  variable  of  integration  to 

Z  =:  p(l-t) 


we  obtain 


1  ”  £  (p)  Pr<1-P)n’r  +  (N-1 ) r (n-N)'!'  / 

r&O  o 

in  which  one  recognises  the  integral  representation  of  the  incomplete  Beta 
function  in  the  last  term,  namely 


(9) 


Bp(a,p) 


But  from  equation  (7)  we  have 


1  -  Pj(n;p) 


Z^O-Z)*3-1  dZ 


N-1 


n\  r/m  Nn-r 


p  0-p) 


r=0 


and  thus,  by  comparison  with  equations  (9)  and  (lO) 


pN(n;p)  =  (N-1)! '(n-N)!  Bp^N>  n_W+1 ) 

B  (N,  n-N+1 ) 

=  T?rT' '  7TTT  =  I  (N,  n-N+1 ) 
B(N,  n-N+1 )  pv  * 


(10) 


(11) 


where  B(a,|3) 
a  and  (3  as 


B,(a,p)  is  the  usual  complete  Beta  function,  defined  for  integer 


)■  • 


(12) 
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The  quantity  I^(a,  (3)  has  been  tabulated  by  Pearson  for  the  range  of 
variables  0  <  p(0*01  )  ^  1  and  a,  |3  <  50.  Using  these  tables  P^(n;p)  has  been 
plotted  in  Figs.1  to  6  for  N  =  3,  6,  9,  12,  15  and  20  respectively  and  for  p 
ranging  from  0*3  to  1*0  in  steps  of  0*1.  For  p  =  1,  of  course,  the  distribution 
is  simply  a  step  function  at  n  =  N. 

These  curves  show  quite  clearly  the  sharp  increase  in  the  number  of 
launches  required  when  the  probability  of  a  successful  launch  is  low  even  when 
no  satellites  fail  between  launches. 

3.3  Gaussian  distribution  approximation  to  P^(n;p) 

It  is  shown  in  Appendix  A  that  for  N  large  the  establishment  distribution 

P®(n;p)  tends  towards  a  Gaussian  distribution  with  a  mean  m  =  N/p  and  variance 
o  2 

CT^  es  N(l-p)/p  •  Thus,  if  it  is  assumed  that  the  distribution  is  Gaussian,  the 
levels  of  probability  for  n  are: 


P(n  $  m  ±43")  = 

P(n  <  m  ±3cr)  = 

P(n  ^  m  ±2cr)  = 

P(n  <  m  ±  cr)  = 

P(n  <  m)  =  50 $ 


50  ±49-9$ 
50  ±49- 85/0 
50  ±47*75$ 
50  ±34*15$ 


(13) 


0  /  \ 

Using  the  Gaussian  distribution  approximation  to  P^(n;pj  the  corresponding 
curves  have  also  been  plotted  for  comparison  with  the  exact  curves  in  Figs.1  to  6, 
for  p  »  0‘5»  From  these  it  is  evident  that  (particularly  for  N  greater  than 
about  10)  the  approximation  is  extremely  good  for  practical  purposes (maximum 
error  for  a  given  probability  is  about  one  launching).  The  main  discrepancy 
between  the  true  and  approximate  distributions  occurs  in  the  50%  probability 
region,  whereas  for  high  degrees  of  probability  the  agreement  is  fairly  good, 
even  for  small  values  of  N.  Consequently,  as  a  useful  working  guide,  it  appears 
that  the  distribution  may  be  adequately  represented  by  (particularly  for  high 
levels  of  probability  )  the  Gaussian  distribution. 

Assuming  a  Gaussian  distribution  approximation,  Fig.  7  shows  the  number  of 
launches  required  to  successfully  establish  N  (<  30)  satellites  in  orbit  with 
probabilities  of  84,  98  and  99*8$  and  with  the  probability  of  a  successful 
launch  ranging  from  p  =  0*3  to  1. 
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MAINTENANCE  OF  THE  SYSTEM  ASSUMING-  BMEDIATE  REPLACEMENT  OF  FAILED 
SATELLITES 


4. 1  Without  launch  failures 

Having  established  a  system  requiring  N  satellites  in  orbit,  let  us  now 
consider  the  maintenance  of  the  system  when  it  is  assumed  that  a  failed 
satellite  may  be  replaced  immediately  by  a  new  one  and,  furthermore,  that  the 
new  launching  is  perfect.  The  effect  of  unsuccessful  launchings  will  be 
considered  later  in  this  section.  Thus,  under  these  ideal  conditions,  the 
number  of  satellites  in  orbit  is  always  maintained  constant  and  we  wish  to 
determine  the  distribution  of  the  firing  rate  demanded  to  achieve  this.  This 
assumption  is  equivalent  to  assuming  ihat  the  time  to  establish  a  satellite  in 
orbit  and  also  the  time  between  possible  successful  firings  are  negligible 
compared  with  satellite's  mean  lifetime. 


Now,  for  a  single  satellite,  since  the  lifetime  has  a  Poisson  distribution, 
the  probability  that  the  satellite  will  be  replaced  exactly  k  times  (immediate 
replacement)  during  a  period  T  is 

p(k?a)  =  fr  e**a  (14) 


where  a  =  T/-6.  This  is  a  well  known  formula  representing,  for  example,  the 

2 

number  of  changes  of  state  for  the  ’telegraph  signal*  «  But  the  sum  of  two 
independent  quantities  having  the  Poisson  distributions  p(k^;a>|)  and  pCk^jct^) 
also  has  the  Poisson  distribution 


pCkj+kgj  Oj+o^)  . 


Consequently,  the  probability  that  n  replacements  will  be  necessary  to  maintain 
N  satellites  in  orbit  for  a  period  T  is 

p(n;oN)  =  e*^  (15) 

which,  of  course,  is  the  same  distribution  as  that  required  for  maintaining  a 
single  satellite  in  orbit  for  a  period  NT#  The  probability  distribution  for 
the  maintenance  of  N  satellites  in  orbit  with  n  launchings  during  a  period  T, 
assuming  immediate  successful  replacement,  is  therefore 


n 

^  p(r;aN) 
r=0 


(4N)r  e-aN 
r! 


r=0 


(oa=T/£)  .  (16) 
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4. 2  Allowing  for  launch  failures 

Let  us  now  assume  that  the  launchings  only  have  a  probability  p  of  being 
successful  but  that  the  failed  satellites  may  still  be  immediately  replaced. 

Thus  the  assumptions  of  para  4*1  still  hold  apart  from  the  proviso  that  several 
instantaneous  launchings  may  have  to  be  made  before  a  successful  replacement  is 
accomplished.  As  before,  the  probability  that  k  satellites  must  be  replaced  is 
p(k;a,N)  given  by  equation  (15).  But  the  probability  that  k  satellites  will  be 
successfully  launched  into  orbit  when  not  more  than  n  are  launched  with  success 
probability  p  is  P^(n;p)  given  by  equation  (6)  or  (7)»  Consequently,  the 
probability  distribution  for  maintaining  N  satellites  in  orbit  with  n  launchings 
over  a  period  T  when  it  is  assumed  that  there  is  immediate  replacement  with 
successful  launching  probability  p  is  simply 

(17) 

or,  explicitly,  using  equations  (7)  and  (15) 

n 

Pjj(n;p)  * 

k=0 

Alternatively,  equation  (6)  may  be  used  instead  of  equation  (7)  giving 

n  n  . 

l\  k,  ,r-k  (oN)  -aN 

ij  p  0-p)  V"  6 

k=0  r=k 


n 


r/„  >.n-r 

P  (1-p) 


k! 


-oN 


(18) 


r=k 


n 


#»»> 


P®(n;p)  p(k;a») 


k=0 


the  individual  terms  of  this  series  now  being  independent  of  n. 

4.3  Evaluation  of  the  probability  distribution  P^(n;p) 

As  we  have  seen  in  section  4*1,  for  p  =  1,  the  distribution  P^  is  determined 

In  x  —a  ' 

by  the  sum  of  a  series  of  Poisson  terms.  Several  tables  of  the  function  — — 
exist  of  which  the  most  comprehensive  appear  to  be  those  compiled  by  Molina^  of 
Bell  Telephone  Laboratories;  these  tabulate  the  function  for  a  <  100  besides 
giving  values  of  the  sum 
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P(c,a)  — 


x  -a 
a  e 


;  P(0,a)  = 


Thus,  comparing  equations  (16)  and  (20),  we  have 


P“(n;l)  =  1  -  P(n+1,  NT/*)  (21 ) 

and  consequently  there  is  no  difficulty  in  plotting  the  distribution. 

The  function  P(c,a)  may  also  be  written  in  terms  of  the  incomplete  Gamma 
function,  since,  if  we  put  f(x)  s  eX,  q  =  -a  and  p  =  a  in  equation  (8),  it  follows 


n-i 


Al-t)1'-1  e-t1-*) 


(N— 1 ) ! 


r  r  (n) 

7T  6  +  rfi^T 


where  r(N)  is  the  complete  Gamma  function  and  T  (N)  is  the  incomplete  Gamma 

a 

function  defined  by 


ra(N) 


I  (l-t)”'1  dt 


a 

-/ 


-t  N-1  ... 
e  t  dt  . 


Consequently,  from  equations  (21 )  and  (22),  it  follows  that  the  desired 
relationship  is 


P(N,a) 


ra(N) 

=  w 


(24) 


in  the  notation  of  Pearson's^  tables  of  the  incomplete  G-amma  function  and  hence, 
using  equations  (21 )  and  (24) 


pJJUji) 


(25) 


For  p  *  1 ,  the  distribution  is  more  difficult  to  calculate  since,  as  is 
seen  by  inspection  of  equations  (18)  or  (19),  it  now  consists  of  a  coupling 
between  the  binomial  and  Poisson  distributions*  Although  P^(njp)  may  be  expressed 
in  terns  of  tabulated  quantities  it  is  an  extremely  tedious  process  to  carry  out 
the  calculations  involved.  For  instance,  P™(n;p)  can  be  calculated  from  the 
expression 


xp(k^  n~k+i ) 

k=0 


-aN 

e 


y 


(26) 


obtained  by  substituting  equation  (ll)  in  (l7)« 

Rather  than  perform  this  calculation  for  many  values  of  n  and  p  a  sim¬ 
plification  was  attempted.  Inverting  the  order  of  summation  in  equation  (18) 
we  have 


n  r 


r=0  k=0 


k! 


-aN 

e 


B  1 


r=0  k=r+1 


n 


1 


n 


pr(l-p)n  r  P(r+1,  aN) 


(27) 
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But,  for  n  large,  the  binomial  distribution  B^(n,p)  becomes  Gaussian^  with  mean 
np  and  variance  np(l-p)#  Consequently,  equation  (27)  may  be  approximated  by 

n 

Pj}(n;p)  *  1  -  ^  GrP(r+1,  oN)  (28) 

r=0 


where  G^  is  the  Gaussian  probability  density 


& 

r 


np(l-p) 


exp 


r  - 

2np(l-p) 


J 


(29) 


In  this  form,  equation  (28)  takes  into  account  the  'tail'  of  the  Poisson 
distribution  for  large  values  of  n.  The  quantity  P(r,  aN)  is  unity  for  r  =  0 
(equation  (20)),  approximately  ^  for  r  equal  to  the  mean  value  aN  of  the 
distribution  and  tends  monotonically  to  zero  for  r  large.  Consequently,  for 
large  values  of  n  the  Gaussian  distribution  G^,  with  mean  np,  lies  on  the  ‘tail’ 
of  P(r,  aN)  and  hence  one  would  expect  equation  (28)  to  give  a  fairly  good 
representation  of  the  distribution  F^(n;p)  for  n  large. 

Figs. 8,  9  and  10  for  N  =  6,  9  and  12  were  constructed  using  this 
approximation  for  F^(  njp).  The  curves  correspond  to  a  range  of  values  of  p 
from  0*5  to  1 *0  and  the  ratio  of  the  system  lifetime  T  to  the  mean  lifetime  & 
of  each  satellite  has  been  taken  equal  to  4»  For  instance,  for  a  mean  lifetime 
per  satellite  of  5  years  the  curves  indicate  the  probability  distribution  of 
the  number  of  launches  required  over  a  period  of  20  years  in  order  to  maintain 
the  system.  For  p  =  1  the  curves  were  obtained  from  tabulated  values  of  P(c,a) 
directly.  Calculations  involving  the  exact  formula  for  P^(n;p)  at  isolated 
points  gave  extremely  good  agreement  with  the  values  obtained  using  the 
approximate  formula. 

4.4  Gaussian  distribution  approximation  to  F^(n,p) 

It  is  well  knownv  that  the  Poisson  frequency  distribution 


( j  =  0,1 ,2  ... 

tends  to  the  Gaussian  distribution  for  large  values  of  aN.  Consequently,  one 
would  expect  that  the  distribution  of  P^,  being  a  combination  of  binomial  and 
Poisson  frequency  distributions  should  also  become  Gaussian  for  aN  large.  That 
this  is  the  case  is  proved  in  Appendix  A  where  it  is  shown  that  the  limiting 
Gaussian  distribution  has  a  mean  value 


and  a  variance  given  by 


cr2  =  2®L£sl  ,  (31) 

P 

The  distribution  curves  for  p  =  0*5  when  it  is  assumed  that  is  normal  with  the 
above  mean  and  variance  are  shown  in  Figs. 8,  9  and  10  for  comparison  with  Ihe 
true  distribution.  It  is  evident  from  these  that  the  normality  assumption  is 
quite  valid  for  practical  purposes  where  one  really  only  requires  an  estimate  of 
the  distribution  curve.  As  expected  the  agreement  is  better  for  large  values  of 
dN.  Fig. 11  also  shows  that,  even  for  the  low  value  of  aN  =  4  and  p  =  1,  the 
discrepancy  between  the  Gaussian  approximation  and  the  exact  distribution  is 
generally  small  and  for  high  levels  of  probability,  ~  1 ,  extremely  good 
agreement  is  obtained. 

Accepting  the  validity  of  the  Gaussian  approximation  the  number  of  satellites 
that  may  be  required  in  order  to  maintain  a  system  over  a  period  of  time  with  a 
given  degree  of  reliability  may  easily  be  constructed.  A  typical  example  of  this 
is  shown  in  Fig.  12  for  T/-6  =  l+f  for  systems  requiring  up  to  30  satellites  in 
orbit  and  p  ranging  from  0*3  to  1 ;  the  curves  have  been  plotted  for  1,  2  and  3  cr 
probability  levels  of  being  able  to  successfully  maintain  the  system. 

5  EFFECT  OF  FINITE  LAUNCH  INTERVAL 

One  of  the  shortcomings  of  the  analysis  carried  out  so  far  is  the  assumption 
that  any  failed  satellite  may  be  replaced  immediately  in  the  maintenance  phase  and 
that  during  the  establishment  phase  the  launchings  take  place  at  infinite  rate 
until  the  system  has  its  correct  complement  of  satellites.  Although  these 
simplifying  assumptions  are  acceptable  under  certain  conditions  it  is  of  interest 
to  know  the  consequences  when  they  are  not  even  approximately  true;  for  example, 
when  there  are  only  a  small  number  of  launch  pads  and  the  time  between  firings 
has  to  be  fairly  large. 

In  order  to  estimate  this  effect  another  model  of  the  system  has  been 
considered,  differing  from  the  previous  model  in  that  launchings  may  only  be 
made  at  intervals  t  both  during  the  establishment  phase  and  the  maintenance  phase. 
Thus  during  the  establishment  phase  an  attempted  launch  is  made  at  every  launch 
interval  until  the  required  number  of  satellites  are  in  orbit.  During  the 
maintenance  phase,  if  the  number  of  satellites  in  orbit  is  less  than  the  required 
number  at  the  time  of  a  possible  launch  then  a  launch  is  attempted;  otherwise  no 
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action  is  taken  until  the  next  possible  launch  time  when  a  similar  decision  to 
launch  is  made.  Since  the  lifetime  of  the  satellites  has  been  assumed  to  have 
a  Poisson  distribution  this  means  that  the  probable  state  of  the  system  at  any 
firing  instant  depends  only  upon  the  state  of  the  system  at  the  previous  possible 
firing  instant.  Prom  a  mathematical  point  of  view  this  means  that  the  system  may 
be  represented  by  a  Markov  chain  with  a  transition  matrix  P  depending  upon  the 
possible  changes  of  state  of  the  system  during  the  launch  interval  v.  Appendix  B 
gives  the  relevant  properties  of  Markov  chains  which  are  required  in  the  following 
analysis.  If  the  ith  state  of  the  system  is  defined  as  that  corresponding  to  i 
satellites  in  orbit  immediately  after  a  firing  then  the  transition  matrix  P  yd.ll  be 
of  order  (N+l),  assuming  the  maximum  number  of  satellites  in  orbit  to  be  N,  with 

elements  p.  .  corresponding  to  the  state  transition  i  to  j. 

J 

In  general  there  are  two  ways  in  which  the  transition  from  state  i  to  j  may 
occur  during  the  time  interval  (i)  either  there  can  be  a  failure  of  (i~ j) 
satellites  and  an  unsuccessful  launch,  or  (ii)  a  failure  of  i-j+1  satellites 
followed  by  a  successful  satellite  launching.  For  cases  in  which  j  is  greater 
than  i+1 ,  the  transition  probability  is  zero,  and  for  those  cases  in  which  j  is 
equal  to  i+1  only  the  transition  of  type  (ii)  is  possible.  Also,  if  there  are  N 
satellites  in  orbit  at  the  time  of  a  possible  firing,  no  new  launch  is  attempted* 
Consequently,  if  f  »  1  -  e  is  the  probability  that  a  single  satellite  will 
fail  during  the  interval  tt,  the  elements  of  the  transition  probability  matrix 
are:- 
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also  be  factorized  into  the 

product  of  two  matrices  in 

the 

P  =  DB 

(34) 

elements  of 

the  matrices  D  and  B  are 

D.  . 
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=  0 
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In  computing^  P  may  conveniently  be  obtained  by  this  matrix  multiplication 
since  the  elements  of  the  ith  row  of  D  are  the  terms  corresponding  to  the 
binomial  expansion  of  (f+s)1  and  consequently  the  elements  of  successive  rows  of 
D  may  be  obtained  by  multiplying  the  previous  row  by  (f+s)#  Also  successive  rows 
of  the  matrix  B  may  be  obtained  by  single  shifts  of  the  elements  of  the  first  row 

( P*  •  •  •  •  o)* 

It  may  also  be  noted  in  passing  that  the  factorization  P  »  DB  gives  a 
matrix  D  depending  only  upon  the  probability  of  failure  f  and  thus  corresponds 
to  the  ’death1  process  of  the  system  whereas  B  only  depends  upon  the  probability 
of  replacement  p  and  thus  represents  the  ’birth*  process  of  the  system  in 
statistical  language* 
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Knowing  the  basic  transition  matrix,  the  probabilities  of1  having  various 
numbers  of  satellites  in  orbit  as  a  function  of  time  may  easily  be  obtained  by 
repeated  multiplication  of  P  by  itself,  since  the  probability  of  having  j 
satellites  in  orbit  after  time  nx  starting  with  i  in  orbit  is  given  by 

P?n)  =  [Pn]..  .  (37) 

ij  u 

A  typical  example  of  these  probabilities  is  shown  in  Fig. 13  which  has  been 
computed  for  the  case  N  =  12,  p  =  0*5  and  f  =  0*02  i. e.  for  a  mean  satellite 
lifetime  of  5  years  this  would  correspond  to  a  launching  rate  of  10/year  since 
f  %/ij  for  f  «  1.  For  purposes  of  comparison,  the  probabilities  of  establish¬ 
ment  for  infinite  firing  rates  (f  =  0)  are  also  shown  in  the  same  diagram  for 
p  =  0*5  and  0*7 .  Quite  clearly,  the  firing  rate  has  a  significant  effect  on 
the  establishment  time.  For  instance,  in  the  case  of  p  =  0*5,  the  finite  firing 
rate  of  10/year  compared  with  the  infinite  rate  increases  by  5Q&  the  number  of 
launches  required  to  establish  the  system  with  95^  probability. 

Since  the  computation  required  by  equation  (37)  is  extremely  time 
consuming  on  a  digital  computer  and  does  not  give  any  detailed  statistics  of 
the  system  an  alternative  approach  was  employed.  We  have  already  seen  that 
the  probability  distributions  appear  to  be  well  represented  by  a  Gaussian 
distribution  and  consequently,  if  we  make  the  assumption  that  the  distribution 
is  approximately  Gaussian  then  the  determination  of  the  mean  and  variance  of  the 
Markov  process  should  give  a  good  indication  of  the  distribution.  Fortunately, 
the  mean  and  variance,  besides  certain  other  statistics,  may  be  obtained  by 
fairly  rapid  procedures  from  the  transition  matrix.  A  summary  of  the  relevant 
formulae  is  given  in  Appendix  B. 

In  order  to  determine  the  mean  and  variance  of  the  number  of  launchings 
for  the  establishment  of  the  system,  one  approach  is  to  make  the  state  N  into 
an  absorbing  state  by  putting  the  elements  p„.  of  the  transition  matrix  equal 
to  zero  for  i  I  N  and  p^  =  1  i.e.  we  ignore  what  happens  to  "the  system  once  it 
has  been  established.  This  modified  transition  matrix  then  generates  an 
absorbing  Markov  chain  for  which  7/e  require  the  mean  and  variance  of  the  number 
of  steps  from  state  0  to  reach  the  absorbing  state  N.  Equations  (b6)  and  (B7) 
of  Appendix  B  provide  this  relevant  information  and  have  been  programmed  for  a 
digital  computer  using  the  fundamental  matrix  N  defined  by  equation  (B2). 

An  alternative  method,  which  also  provides  additional  useful  information, 
is  to  use  the  formulae  given  in  Appendix  B  for  the  regular  Markov  chain 
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generated  by  the  original  transition  matrix  P.  Here  the  fundamental  matrix  Z 
is  defined  by  equation  (B9)  and  is  used  in  equation  (BIO)  to  give  the  matrix  M 
with  elements  m^.,nn  ^  being  the  mean  number  of  firings  starting  in  state  i 
before  reaching  the  state  j  for  the  first  time.  In  particular  m^,  is  'the  mean 
number  of  firings  required  to  establish  j  satellites  and  clearly 


m.  . 


=  m.,  + 

xk 


m. 


kj 


for  k  lying  between  i  and  j.  The  variance  of  the  number  of  firings  for  establish¬ 
ment  is  given  by  equation  (B12).  It  is  also  worth  noting  that  if  (N+1 )  is  the 
order  of  the  matrix  P  then  the  mean  and  variance  of  the  required  number  of  firings 
for  given  f  and  p,  is  obtained  for  the  establishment  of  all  systems  requiring  no 
more  than  N  satellites  in  orbit.  Finally,  the  elements  (aQ,  a^  ....  a^)  of  the 
row  vector  a  which  is  the  solution  of  aP  =  a,  represent  the  fraction  of  the  time 
that  the  satellite  system  can  be  expected  to  be  in  the  state  of  0,  1,  ...  or  N 
satellites  in  orbit  over  a  long  period  of  time.  Consequently,  this  provides 
information  regarding  the  likely  state  of  the  system  over  the  maintenance  phase. 

A  computer  programme  has  been  written  to  calculate  these  various 
quantities.  Although  only  the  case  of  single  launchings  has  so  far  been 
investigated  the  programme  was  written  for  the  more  general  case  of  multiple 
launchings  by  appropriate  modification  to  the  transition  matrix  P. 

Figs. 15,  16  and  17  show  the  expected  number  of  launches  which  would  have 
to  be  made  to  establish  systems  of  6,  12  and  20  satellites  with  probabilities 
50,  84,  98  and  99*  8/S  as  a  function  of  the  launch  probability  p  and  the  launch 
interval  parameter  f  ^  t/£,  assuming  the  Gaussian  distribution,  (For  reference 
purposes  the  quantity  f  has  been  plotted  in  Fig.  1 4.  as  a  function  of  the  launch 
interval  for  various  mean  satellite  lifetimes. )  It  is  clear  from  these  diagrams 
that  the  launch  rate  must  be  kept  quite  high,  even  for  systems  requiring  a  few 
satellites,  if  the  probability  of  successfully  launching  a  satellite  is  not 
high,  otherwise  a  very  large  number  of  launchings  would  be  required.  This  is 
simply  due  to  the  fact  that  there  is  a  high  chance  of  a  satellite  failing 
during  the  mean  interval  i/p  between  successful  launches.  For  instance,  for  a 
system  consisting  of  6  satellites  with  mean  lifetimes  of  5  years,  if  the  launch 
interval  is  3  months  i.e.  f  =  OO5  and  the  probability  of  a  successful  launch 
is  0*5,  the  expected  time  between  successful  launchings  is  3/0*5  =  6  months# 
Consequently  the  probability  of  at  least  one  of  the  satellites  failing  during 
this  time  is  Pf  =  1  -  (l-f')6  where  f’  =  f/0*5  =  0*1  so  that  P  =  4 jfo  and  it 
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is  not  surprising  that  40  launches  i.e.  10  years  (Fig.15)  are  needed  to 
establish  the  system  with  99*8$  confidence.  On  the  other  hand,  if  the 
firing  interval  is  halved,  i.  e.  1-$  months  and  f  =  0*025,  30  launchings  over 
nearly  4  years  are  now  required. 

The  only  real  remedy  for  establishing  a  system,  for  a  given  firing  rate,  in 
a  reasonably  short  time  is  to  increase  the  reliability  of  the  launch  or  the 
mean  lifetime  cf  the  satellites.  For  instance,  if  p  =  0*7  in  the  last  example 
above  then  the  number  of  launches  is  reduced  from  30  to  about  18  over  2  years. 

But  if  the  launch  rate  is  high,  i.e.  f  small,  there  is  little  to  be  gained,  at 
least  for  the  establishment  phase,  by  increasing  the  mean  lifetime,  as  is 
indicated  by  the  steep  slope  of  the  curves  in  Figs. 15,  16  and  17»  This  is  as 
expected  since  there  would  now  be  very  few  satellite  failures  during  the 
establishment  period.  In  fact,  although  the  number  of  launches  must  be  greater 
for  f  *  0,  the  effect,  as  may  be  seen  from  the  curves  in  Fig.15,  16  and  17, 
only  becomes  really  significant  as  f  is  increased  above  a  certain  value.  On 
the  other  hand,  the  main  disadvantage  of  increasing  r ,  whilst  retaining  f 
constant  by  increasing  the  satellite  lifetime,  is  that  the  establishment  time 
is  increased  and  this  is  usually  of  prime  importance. 

For  instance,  let  us  assume  that  a  system  of  12  satellites  is  required 
in  orbit,  the  mean  lifetime  of  each  satellite  being  5  years.  If  there  are  4 
launch  pads,  each  capable  of  firing  a  satellite  every  3  months  with  a  success 
probability  of  0*7(v  =  l/l6  year,  f  =  0*0125)  then  (Fig.l6)  the  expected  number 
of  firings  would  be  about  26  over  a  period  of  If  years  in  order  to  establish  the 
system  with  a  probability  of  98%.  Increasing  the  lifetime  to  10  years 
(f  =  0*00625)  would  save  only  about  2  launches  although,  if  the  mean  satellite 
lifetime  was  only  2^  years  (f  =  0*025),  an  increase  of  8  firings  would  be 
expected.  On  the  other  hand,  if  the  firing  rate  was  halved  (^  =  5  years, 
f  =  0*025)  e. g.  2  pads  each  launching  at  4  per  year,  the  establishment  time 
would  be  increased  to  4?  years!  Lower  launch  success  probabilities  would  make 
these  changes  even  more  dramatic. 

Clearly,  very  little  can  be  said  here  regarding  the  best  manner  in  which 
to  improve  the  system  in  general  since  it  depends  upon  the  state  of  the  various 
parameters  •&,  p,  f,  N  and  the  cost  of  changing  any  of  these  quantities.  Even 
so  the  various  curves  allow  one  to  determine  the  time  of  establishment  of  a 
system,  given  the  various  parameters  for  a  particular  launch  scheme  and  thus 
to  judge  whether  such  a  process  is  feasible. 
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After  establishment.  Figs. 18  to  22  show,  as  a  function  of  f,  the  percentage 
of  time  that  there  are  i  satellites  in  orbit  for  systems  of  (nominally)  6,  12 
and  20  satellites.  Although  the  objective  is  to  maintain  N  satellites  continuously 
in  orbit  there  is  always  a  finite  probability  that  this  will  not  be  achieved 
unless  the  firing  rate  is  infinitely  fast  (f  =  0).  For  instance,  for  a  system  of 
6  satellites  and  p  =  0*7  (Fig. 18)  if  f  is  less  than  0*02  there  is  negligible 
probability  that  the  number  of  satellites  in  orbit  would  drop  below  5  and  the 
probability  of  6  satellites  is  greater  than  86%.  Consequently,  for  f  =  0*02, 
it  may  be  concluded  that  the  corresponding  firing  rate  would  be  quite  suitable 
for  maintaining  a  system  of  5  satellites  (although  still  attempting  to  maintain 
6  satellites)  and  similarly  f  ==  0°04  would  maintain  4  satellites  and  so  on.  On 
the  other  hand,  if  p  was  only  0*5,  then  f  =  0*02  would  only  be  suitable  for 
maintaining  4  satellites  and  f  =  0*04  for  3  satellites.  Conversely,  for  a  given 
f,  one  may  determine  the  number  of  satellites  which  one  should  attempt  to 
maintain  in  orbit  so  that  at  least  a  certain  number  should  be  maintained  with  a 
given  probability. 

For  example,  if  f  =  0*01  and  p  =  0«7,  then  from  Fig. 22(a)  in  order  to 
maintain  18  or  more  satellites  in  orbit  a  possible  strategy  would  be  to  attempt 
to  retain  20  satellites  functioning.  The  20  satellites  would  be  maintained  with 
a  probability  of  ~!0>,  at  least  19  with  94%  and  at  least  18  with  a  probability  of 
98%  and  negligible  probability  of  these  being  less  than  18. 

Finally  it  is  worth  mentioning  two  facts  which  have  not  been  considered  in 
this  Report  and  which  may  have  a  significant  effect.  In  the  first  place  no 
account  has  been  taken  of  the  time  taken  to  place  a  satellite  in  orbit  from  its 
time  of  firing.  Secondly,  the  fact  has  not  been  included  that  a  launching  may 
be  attempted  at  any  time,  rather  than  at  discrete  intervals,  after  the  launch 
pad  preparation  time  has  elapsed  from  its  last  firing.  In  order  to  assess  these 
effects  an  entirely  new  theoretical  approach  is  required,  analogous  to  the  theory 
of  the  maintenance  of  machines,  but  so  far  a  satisfactory  theory  has  not  been 
found. 

6  CONCLUSION 

During  the  establishment  phase  of  a  satellite  system,  if  the  number  of 
satellite  failures  is  negligible,  i.e.  the  launching  interval  is  very  small 
compared  with  the  satellite  mean  lifetime  (t/4  ^  0),  then  the  probability 
distribution  of  establishing  the  system  of  N  satellites  is  approximately 

2 

Gaussian  (for  practical  purposes)  with  a  mean  of  N/p  and  variance  N(l-p)/p  . 
Consequently  a  fairly  rapid  estimate  of  the  likely  number  of  launches  for  a 
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given  probability  of  establishment  may  be  obtained  without  evaluating  the  exact 
distributions  which  have  been  given  in  this  Report.  As  the  launch  success 
probability  is  decreased  there  is  a  drastic  increase  in  the  number  of  required 
launches  due  to  the  mean  of  the  distribution  being  inversely  proportional  to  p. 

The  exact  nature  of  this  increase  is  illustrated  in  Figs.1  to  7  as  a  function 
of  the  number  of  satellites  required  in  orbit. 

A  finite  time  interval  between  launches  only  affects  the  number  of  launches 

-v/A 

for  establishment  significantly  when  the  ratio  t / A,  and  hence  f  =  1-e  ,  is 

increased  above  a  certain  region.  This  is  indicated  by  the  steep  slope  of  the 
curves  near  f  =  0  in  Figs. 15  to  1 7  and  is  due  to  the  increased  probability  of 
satellite  failures  as  the  time  interval  between  launches  is  increased  or  the 
mean  lifetime  is  decreased.  For  a  particular  satellite  system  the  effect  of  varying 
the  launch  interval,  probability  of  a  successful  launch  or  mean  lifetime  on  the 
number  of  launches  or  the  establishment  time  may  be  obtained  from  the  curves  in 
these  figures. 

During  the  maintenance  phase,  if  it  is  assumed  that  any  failed  satellites 
may  be  replaced  immediately  (f  =  0),  then  the  probability  distribution  of  the 
number  of  firings  over  a  period  T  is  also  approximately  Gaussian  with  a  mean  of 

O 

aN/p  and  variance  aN(2-p)/p  where  a  =  T/-6.  Using  this  distribution  simple 
calculations  may  be  made  to  estimate  the  required  firing  rate  after  establish¬ 
ment.  Low  values  of  p  will,  of  course,  have  a  marked  effect  on  the  launching 
rate.  This  is  illustrated  in  Fig.1 2  where  the  number  of  launches  for  maintain¬ 
ing  a  system  of  N  satellites  (N  <  30  )  with  a  given  level  of  probability  and  a  =  4 
is  shown  for  various  values  of  launch  success  probability. 

Using  a  launch  strategy  which  attempts  to  maintain  a  specific  number  of 
satellites  in  orbit  so  that  launches  are  only  made  when  the  number  of  satellites 
in  orbit  falls  below  this  number,  gives  rise  to  the  probability  distribution  of 
the  percentage  of  time  that  there  are  more  than  a  certain  number  of  satellites 
in  orbit.  Consequently,  by  attempting  to  maintain  a  number  of  satellites  in 
orbit  which  is  in  excess  of  the  number  required  for  communication  purposes  if 
there  were  no  failures,  the  probability  of  breaks  in  communication  occurring  in 
the  actual  system  may  be  decreased.  Estimates  of  the  surplus  number  of  satellites 
which  must  be  held  in  orbit  so  that  a  system  may  function  as  a  communication 
system  with  little  loss  in  transmission  time  have  been  discussed  in  section  5« 

One  or  two  excess  satellites  appear  to  be  sufficient,  for  satellite  systems 
requiring  less  than  20  satellites  provided  the  launch  success  probability  is 
greater  than  7Q&  and  f  less  than  0*01. 
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MOMENTS  AM)  LIMITING  DISTRIBUTIONS 


The  moment  generating  function  M(s)  of  the  establishment  distribution, 
equation  (5),  is  given  by 


OO 


r=0 


where  Li*  is  the  rth  moment  about  the  origin  of  the  distribution  pfT(n.p).  The 
rr  0  Nx  r 

moments  about  the  mean,  (i  ,  are  obtained  from  [a*  by  means  of  the  relationship 


or 


K 


Z_7  \ 

j=0 


r)  |i'  .  (-H»  )' 

J/  rr~j 


j=0 


(A3) 


(A4) 


In  particular, 


=  m2 


=  M3  "  3  H2  +  2  n1 


‘3 


(A3) 

(AO 


and 
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It  also  follows  by  differentiation  of  equation  (A2)  that 


t 


[£, 

-ds1 


» 


s=0 


<A7) 


Thus,  from  equation  (A1 ) , 


M(s)  = 


-s 

-e  -  q~ 


N 


M(0)  =  1 


and  therefore 


M1  (  s 


Ne 


-s 


N 


-s  .  s 

e  -  q  1-qe 


(AS) 


Hence,  the  mean  is 


=  M«(0)  =  |  . 


(A9) 


Differentiating  (A8), 


(1-  qeS)  M”  (s)  -  qM«  (s)  =  Mlf(s) 


i.e, 


(1  -qeS)  M“  (s)  =  (N+q)  M*(s) 


(A10) 


thus  gives  the  second  order  moment  about  the  origin  as 


(i*,  =  M"  (0)  = 


2 
P 


and  hence  the  variance  of  the  distribution  is 


O'  =  H, 


_u>2  _  rc(Mfi-p) 
2 


If 


i»e4 


^  -  2^  . 
P 


(All) 
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Knowing  the  mean  and  variance,  the  moment  generating  function  of  the 

distribution  in  standard  measure  M  (s)  may  be  obtained  by  multiplying  M(s)  by 

11 

***  o-  5 

e  after  replacing  s  by  s/o  *  Thus,  using  equations  (A9)  and  (All) 


.  siNq 
1  -qe  \ 


N 


6  N 


1 


-N 


and  the  cumulative  function 

K  (s)  =  log  M  (s) 

n  D  n 


2 

3 

Thus  for  any  finite  s,  the  cumulative  function  tends  to  «— ■  as  N  •»  <x>,  i*e*  the 
distribution  becomes  Gaussian  for  N  large  since  all  the  cumulants  K  for 
r  >  2  tend  to  zero* 

A* 2  Mai.ntenance  phase 

The  probability  of  maintaining  the  system  with  exactly  n  launchings 
(in  the  time  T  =  a  &  )  is 
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p“(n;p)  =  P”(n;p)  .P”(n-1.p) 

and  thus  with  the  aid  of  equation  (19)*  which  is  more  convenient  than 
equation  (18)  since  only  the  summations  depend  upon  the  variable  n, 

(A12) 

Consequently,  the  corresponding  moment  generating  function  is 

oo 

M(s)  =  ^  pJJ(n;p)  eSn 
n=0 


co  n 


n=0  k=0 


pk(l-p) 


n-k 


(aN)k  -aN  sn 
yf—  e  e 


oo 

V 

L  k! 

k=0 


hxN 

e 


I'pfru-p).’  ~i 

v-py  Li  _  (i-p)esJ 


k 


-ctN 

e 


aNpe 

1  -  (l-p)eS 


x*e, 


»  .  exp  f- 

t  1  -  0-p)< 


(A13) 


Thus 


(l  -qeS)  log  M(s)  =  -  aN(l-eS) 


and  differentiating. 


(1  -  qeS)  -  qeS  log  ll(s)  =  aNeS  * 


(A1A) 
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On  putting  s  =  0  in  equation  (A14)  and  noting  that  M(o)  =  1,  then  the  mean  of 
the  distribution  is 

H*  «  H'(o)  =  —  .  (A15) 

Rewriting  equation  (A14)  in  the  form 


(e  S  -  q)  Jll(s)  -  q  M(s)  log  M(s)  =  aN  M(s) 


and  differentiating  once  again  gives 


(e  S  -  q)  M"(s)  -  e  S  M*  (s)  -  q  M* (s)  log  M(s)  -  q  MT(s)  =  aN  Ml(s)  . 

Consequently,  on  putting  s  =  0,  and  using  equation  (A15)»  we  have 

p  M"(o)  a  — T  (aN  +  1  +  q)  • 

P 

Henoe 

[j.*  =  M"(0)  =  ~  (aN  +  2  -  p)  (A16) 

P 

and  the  varianoe  of  the  distribution  is 


cr 


2 


P 


(A17) 


As  in  the  previous  section,  the  moment  generating  function  of  the 
distribution  in  standard  measure  may  now  be  found  using  the  mean  and  varianoe 
as  given  by  equations  (A15)  and  (A17)»  Thus  we  have 


and,  consequently,  the  cumulative  function  is 


Appendix  A 


27 


and  the  distribution  becomes  Gaussian  as  N  *►  oo. 

It  is  also  to  be  noted  that,  since  a  in  equation  (13)  only  appears  in  the 
form  aN,  the  distribution  will  tend  to  become  Gaussian  as  a  •»  «  also. 
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STATISTICS  OF  A  HARKOV  CHAIN 


B.1  Basic  concepts 


A  Markov  chain  is  defined  in  terms  of  the  transition  probabilities  p^(n) 
between  state  s^  at  step  n  and  state  s,.  at  step  (n+1 )  where  P.y(n)  is  independent 
of  n  and  the  outoomes  before  step  n,  and  may  therefore  be  denoted  by  p. .. 
Consequently,  the  probability  of  being  in  state  s  .  after  n  steps  starting  in 
state  s^  is  simply 


,(n) 


p:.'  =  [p*|.. 

iJ  10 


(B1) 


where  P  is  the  transition  matrix  with  elements  p.  . . 

10 

The  states  of  a  Markov  chain  may  be  divided  into  transient  and  ergodic  sets# 
The  former  sets,  once  left,  are  never  entered  again,  while  the  latter,  once 
entered,  are  never  again  left.  Furthermore,  if  a  state  s^  is  the  only  element  of 
an  ergodic  set,  then  it  is  called  an  absorbing  state  and  will  have  the  element 
p..  =1  with  all  other  entries  of  this  row  p.  .  =0#  If  the  only  non- transient 
states  of  a  chain  are  absorbing  it  is  called  an  absorbing  chain.  On  the  other 
hand,  a  regular  Markov  ohain  is  such  that  it  is  possible  to  reach  any  state 
after  n  step  regardless  of  the  initial  state  and  thus  contains  a  single  ergodio 
set  and  no  transient  sets. 

The  matrix  P  given  by  equation  (33)  in  the  main  text  is  regular  since 
whatever  the  number  of  satellites  in  orbit  at  one  instant  there  is  always  a 
finite  probability  of  getting  into  any  other  state  at  a  later  time.  On  the  other 
hand,  if  one  was  only  interested  in  the  establishment  phase,  the  elements  p.„ 

could  be  put  equal  to  zero  for  i  t  N  and  p^  =  1  so  that  the  state  N  was 

absorbing  and  thus  P  would  generate  an  absorbing  Markov  chain.  This  is  useful 
when  one  is  not  interested  in  what  happens  after  reaching  state  N  i.e.  we  may 

consider  the  attainment  of  state  N  as  the  end  of  the  process. 


B,2  Absorbing  Markov  chains 

In  what  follows,  one  will  only  quote  some  of  the  more  useful  results 
required  in  analysing  Markov  chains.  For  a  proof  of  the  formulae  the  reader  is 

g 

referred  to  Kemeny  and  Snell  , 


Let  (Var.)  denote  the  mean  (variance)  of  a  function  starting  from  the 

state  i  and  n.  the  total  number  of  times  in  the  transient  state  s.*  Also  let 
0  J 
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Q  be  the  transient  part  of  the  transition  matrix  P  obtained  by  deleting  all 
those  rows  and  columns  containing  a  unit  element.  Then,  for  an  absorbing 
Markov  chain,  the  fundamental  matrix  N  is  defined  by 

M  =  (I  -  Q)"1  I 


=  N  . 


iM,  [n|j  =  11(2  -  I) 


where  is  the  matrix  obtained  from  N  by  putting  all  off  diagonal  elements 
equal  to  zero.  Consequently,  the  variance  of  n .  is 


[  Var±[  n j] }  =  N(2Ndg  -  i)  -  N 


where  N  is  the  matrix  obtained  from  N  by  squaring  each  element, 
sq 

Finally,  if  t  is  the  number  of  steps  in  which  the  process  is  in  a 
transient  state,  then 


{M±[t]]  =  Ng  h  v  , 


{ M±[  t  ]  ]  =  (2N-X)t 


so  that 


{Varjt]}*  (2N-i)t-t 


where  g  denotes  a  oolumn  matrix  with  unit  elements. 

For  the  establishment  phase  of  satellites  where  the  Nth  state  is  con¬ 
sidered  absorbing,  equations  (B 6)  and  (B7)  thus  give  the  mean  and  variance  of 
the  time  required  to  establish  the  system. 

B3  Regular  Markov  chains 

If  P  is  regular  then  Pn  tends  to  a  oonstant  probability  matrix  A  as  n  ■+  oo 
with  all  its  rows  equal  to  the  same  probability  vector  a  =  (aQ,a^,  . ...  a^)  and 

a  P  =  a  .  (B8) 
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This  means  that  any  long  range  predictions  are  independent  of  the  initial  state. 
Also  a.  represents  the  fraction  of  time  that  the  system  can  be  expected  to  be  in 
state  s.  over  a  long  period  of  time. 

For  regular  chains,  the  so  called  fundamental  matrix  is  defined  by 

Z  =  [i  -  (P-A)]**1  (B9) 


and  is  the  basic  quantity  used  to  compute  most  of  the  interesting  descriptive 
quantities. 


If  f  .  denotes  the  number  of  steps  before  entering- state  s.  for  the  first 
J  J 

time,  and  if  we  denote  [ll.[f .]]  by  the  matrix  M  =  [m.  .]  ,  M  is  given  by 


M  =  (I  -  2  +  E  Z,  )fi  (BIO) 

dg 

where  R  is  the  diagonal  matrix  with  elements  r. ,  »  i/a. ,  =  m.  ,  and  E  is  a 

IX  IX  XX 

square  matrix  composed  of  unit  entries.  Thus  M  is  a  measure  of  the  mean  number 
of  steps  to  the  first  passage  through  any  state.  In  particular  m  .  is  the  mean 

O  J 

number  of  launchings  to  establish  j  satellites  in  orbit  in  the  context  of  this 
Report. 

The  variance  of  the  first  passage  time  may  also  be  obtained  since,  if 

W  a  f  M.f  f^] }  ,  then  W  is  given  by 
1  J 

W  =  M(2Zdg  R  -  I)  +  2  (ZM  -  E(ZM)dg)  (B1 1 ) 


and  consequently 


Var.[  f .] 
i  J 


sq 


(B12) 


W  -  M 
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a. 

l 

A 

B 

BN(n,p) 

B(a,p) 

Bp(a,p) 

D 

E 

f  =  1-s 

f . 

J 

&r 

I 

l(a,P) 

I  (a,P) 

k!=) 


m 


m.  . 


M 

M. 

1 

M(s) 

n 


n . 
0 
N 


P 


Pj}(n;p) 

p|(n;p) 

p(kja) 


P 

P(t) 


elements  of  vector  a 

matrix  with  all  rows  equal  to  the  limiting  probability  vector  a 

'birth*  matrix,  see  equation  (36) 

binomial  distribution,  see  equation  (4) 

complete  Beta  function  =  B^(a,p) 

incomplete  Beta  function 

'death*  matrix,  see  equation  (35) 

square  matrix  composed  of  unit  entries 

satellite  failure  probability  during  interval  x 

number  of  steps  before  reaching  state  j  for  the  first  time 

Gaussian  distribution,  see  equation  (29) 

unit  matrix 

incomplete  Gamma  function  (Pearson  tables),  see  equation  (23) 
incomplete  Beta  function  (Pearson  tables),  see  equation  (ll) 
cumulative  generating  function 
mean  lifetime  of  a  satellite 
mean  of  a  distribution 

mean  number  of  steps  starting  in  state  i  before  reaching  state  j 
for  the  first  time 

matrix  with  eLements  nn  ^  defined  by  equation  (BIO) 
mean  of  a  function  starting  in  state  i 
moment  generating  function 
number  of  launches 

number  of  times  in  a  transient  state  j 

number  of  satellites  in  complete  system;  fundamental  matrix  of  an 
absorbing  Markov  ohain,  see  equation  (B2) 
probability  of  a  successful  launch 
elements  of  transition  matrix  P 

Poisson  probability  density  of  satellite  lifetime,  see  equation  (l) 

probability  density  for  establishment 

probability  density  for  maintenance 

Poisson  frequency  distribution,  see  equation  (14) 

transition  matrix 

see  equation  (2) 

probability  distribution  for  establishment 
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p>;p) 

P(o,a) 

q 

Q 

R 

s 

t 

T 

Var. 

x 

V/ 

z 

a 


probability  distribution  for  maintenance 
see  equation  (20) 

1“P 

transient  part  of  the  transition  matrix 
diagonal  matrix  Tilth  elements  in. 

1-f,  satellite  survival  probability  during  interval  x 
time 

number  of  launches  before  system  established 
maintenance  period 

variance  of  a  function  starting  in  state  i 
see  equation  (Bll) 

fundamental  matrix  for  a  regular  Markov  ohain,  see  equation  (B9) 
T/£,  see  section  4 

or  limiting  probability  vector  a ,  see  equation  (B8) 

complete  G-amma  function 

incomplete  G-amma  function 

rth  moment  about  the  origin 

rth  moment  about  the  mean 

column  matrix  with  unit  entries 

variance  of  a  distribution 

interval  between  launches 
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Fig.1 


FIG.  I  ESTABLISHMENT  PROBABILITY  FOR  HIGH  LAUNCH  RATES  (N=3) 
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FIG.  4  ESTABLISHMENT  PROBABILITY  FOR  HIGH  LAUNCH  RATES  (N  =  12) 
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FIG.  8  MAINTENANCE  PROBABILITY  DISTR 
FOR  N  =  6  AND  =  4  («N=  24) 
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Fig. 11 


FIG.  1 1  MAINTENANCE  PROBABILITY  DISTRIBUTION 
FOR  N  =  I  AND  T/f  =  4  (a.N  =  4) 


Fig. 12 
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FIG.  13  PROBABILITY  OF  NUMBER  OF  SATELLITES  IN  ORBIT  AFTER  n  LAUNCHINGS 

(n=I2,  p-O-5,  f-0-02') 
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FIG.  14  PROBABILITY  OF  SATELLITE  FAILURE  (/)  AS  A 
FUNCTION  OF  LAUNCH  PERIOD  (t) AND 
SATELLITE  MEAN  LIFETIME  (/>) 
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Fig. 16 
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FIG. 16  NUMBER  OF  LAUNCHES  FOR  THE  ESTABLISHMENT 
OF  12  SATELLITES  AS  A  FUNCTION  OF  / 


Fig. 17 


FIG.  17  NUMBER  OF  LAUNCHI 
OF  20  SATELLITES  A 
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Fig. 18 
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FIG.  18  STATE  DURATIONS  AS  A  FUNCTION 
OF  /  FOR  N  =  6  AND  b  =  0-7 
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FIG.  20  STATE  DURATIONS  AS  A  FUNCTION 
OF  /  FOR  N  =  1 2  AND  f>  =  0-7 
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FI  G.  2 1  -  STATE  DURATIONS  AS  A  FUNCTION 
OF  /  FOR  N  *12  AND  b  -  0-5 
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Fig. 22 
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FIG 22  (cub)  STATE  DURATIONS  AS  A  FUNCTION  OF  /  FOR  N=20,  p--0-7  (a)  AND  p=0-5  (b) 
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